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Abstract: In this talk we discuss how deeply the region of high parton densities has been 
studied experimentally at HERA. We show that the measurements of deep inelastic structure 
functions at HERA confirm our theoretical expectation that at HERA we face a challenging 
problem of understanding a new system of partons: quarks and gluons at short distances with 
so large densities that we cannot treat this system perturbatively. We collect all experimental 
indications and manifestations of specific properties of high parton density QCD. 
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1 What Are Shadowing Corrections? 



In the region of low x and low we face two challenging problems which have to be resolved in 
QCD: 

1. Matching of "hard" processes, that can be successfully described in perturbative QCD 
(pQCD), and "soft" processes, that should be described in non-perturbative QCD (npQCD), 
but actually, we have only a phenomenological approach for them; 

2. Theoretical approach for the high parton density QCD (lidQCD) which we reach in the 
deep inelastic scattering at low x but at sufficiently high Q^. In this kinematic region we 
expect that the typical distances will be small but the parton density will be so large that 
a new non perturbative approach shall be developed for understanding this system. 

We are going to advocate the idea that these two problems are correlated and the system of 
partons always passes the stage of hdQCD before ( at shorter distances ) it goes to the black 
box, which we call non-perturbative QCD and which, practically, we describe in old fashion 
Reggeon phenomenology. In spite of the fact that there are many reasons to believe that such 
a phenomenology could be even correct, the distance between Reggeon approach and QCD is 
so large we are loosing any taste of theory doing this phenomenology. In hdQCD we still have 
a small parameter ( running QCD coupling as ) and we can start to approach this problem 
using the developed methods of pQCD However, we should realize that the kernel of the hd 
QCD problems are non-perturbative one, and therefore, approaching hdQCD theoretically we are 
preparing a new training grounds for searching methods for npQCD. 

First, let me recall that DIS experiment is nothing more than a microscope and we have two 
variables to describe its work. The first one is the resolution of the microscope, namely, Aa; ~ 
1/Q where is the virtuality of the photon. It means that out microscope can see all constituents 
inside a target with the size larger that Ax. The second variable is time of observation. It 
sounds strange that we have this new variable, which we do not use, working with a medical 
microscope. However, we are dealing here with the relativistic system which can produce hadrons 
(partons). So, for everyday analogy, we should consider rather a box with files which multiply 
and their number is, certainly, different in different moment of time. To estimate this time 
we can use the uncertainty principle At oc 1/ AE where AE is the change of energy, namely, 
AE" = Einitiai — E final, and for system of quark and antiquark AE = qo — Qi — 12 = Qo — Q = 
(go-gKgo+g)?" _ q_rn _ where m is mass of the target and go and q is the energy and momentum 

of the virtual photon. Finally, t = 1/ mx with x = ^ where W is the energy of photon - target 
interaction. 

Therefore, the question, that we are asking in DIS at low x, is what happens with constituents of 
rather small size after long time. It is clear that the number of these constituents should increase 
since in QCD each parton can decay in two partons with the probability Pi = i^aas. ^ where 

Ei and ki are energy and momentum of an emitted parton i. 
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This growth we can describe introducing so called structure function ( xG{x, Q^) ) or the number 
of partons that can be resolved with the microscope with definite and x. Indeed, 



xG{x,Q') . (1) 



ain(l/x)91nQ2 7r 

This equation is the DGLAPp] evolution equation in the region of low x. It has an obvious 

solution xG{x, Q^) oc exp {2^ ^"^^ ln(l/x) \ti{Q'^ / Ql)^ . Therefore, we expect the increase of the 
parton densities at a; — > 0. 

In Fig. 1 we picture the parton distributions in the transverse plane. At x ~ 1 there are several 
partons of a small size. The distance between partons is much larger than their size and we can 
neglect interactions between them. 

However, at x — the number of partons becomes so large that they are populated densely 
in the area of a target. In this case, you cannot neglect the interactions between them which 
was omitted in the evolution equations ( see Eq. ( |1] ) ). Therefore, at low x we have a more 
complex problem of taking into account both emission and rescatterings of partons. Since the 
most important in QCD is the three parton interaction, the processes of rescattering is actually 
a process of annihilation in which one parton is created out of two partons (gluons). 

Therefore, at low x we have two processes 

1. Emission induced by the QCD vertex G + G ^ G with probability which is proportional to 
as p where p is the parton density in the transverse plane , namely 

xG^c^ 

where nR^ is the target area; 

2. Annihilation induced by the same vertex G + G ^ G with probability which is proportional 
to asCTQ where as is probability of the processes G + G ^ G, do is the cross section of 
two parton interaction and ctq oc p gives the probability for two partons to meet and 
to interact, while a^Uo p^ gives the probability of the annihilation process. 

Finally, the change of parton density is equal to ^ [H 

9V(x,Q2) N,as -2 



Pi^.Q') - ^p'{x,Q') (3) 



(91n(l/a;)(91nQ2 ^ ry ,^ ; 

or in terms of the gluon structure function 

d'^xG{x,Q'^) Ncas (^s1 ( r^, ^2n\2 



ain(l/x)91nQ2 vr ' ' t^R^Q 

where 7 has been calculated in pQCD | 



xG{xM^) - (xG(a;,Q2)) , (4) 
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Figure 1: Parton distribution in the transverse plane 



3 



Therefore, Eq.( ^ ) is a natural generalization of the DGLAP evolution equations. The question 
arises, why we call shadowing and /or screening corrections such a natural equation for a balance 
of partons due to two competing processes. To understand this let us consider the interaction of 
the fast hadron with the virtual photon at the rest ( Bjorken frame ). In the parton model, only 
the slowest ( "wee" ) partons interact with the photon. If the number of the "wee" partons 
is not large, the cross section is equal to aoN. However, if we have two "wee" partons with the 
same energies and momenta, we overestimate the value of the total cross section using the above 
formula. Indeed, the total cross section counts only the number of interactions and, therefore, 
in the situation when one parton is situated just behind another we do not need to count the 
interaction of the second parton if we have taken into account the interaction of the first one. It 
means that the cross section is equal to 

atot = croN{ 1 - ^ } , (5) 

where R is the hadron radius. One can see that we reproduce Eq.( ^ ) by taking into account 
that there is a probability for a parton not to interact being in a shadow of the second parton. 



2 What Have We Learned about SC? 

During the past two decades high parton density QCD has been under the close investigation of 
many theorists [0] @ @1 0] and we summarize here the result of their activity. 

• The parameter which controls the strength of SC has been found and it is equal to 

- = ^-—^ = -oXp(a:,Q^). (6) 

The meaning of this parameter is very simple. It gives the probability of interaction for two 
partons in the parton cascade or, better to say, a packing factor for partons in the parton 
cascade. 

• We know the correct degrees of freedom at high energies: colour dipoles [^]. By definition, 
the correct degrees of freedom is a set of quantum numbers which mark the wave function 
that is diagonal with respect to interaction matrix. Therefore, we know that the size and 
the energy of the colour dipole are not changed by the high energy QCD interaction. 

• A new scale Q'^{x) for hdQCD has been traced in pQCD approach which is a solution to 
the equation 

Stt^Qs ^ xG{x,Ql{x)) _ 

- mx) Vr^ - ' • 

This new scale leads to the effective Lagrangian approach which gives us a general non- 
perturbative method to study hdCD. 
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• We know that the GLR equation ( see Eqs.( ^ ) and ( ||) ) describes the evolution of the 
dipole density in the full kinematic region We understood that the Mueller-Glauber 
approach for colour dipole rescattering gives the initial condition to the GLR equation. 

• The new, non-perturbative approach, based on the effective Lagrangian P], have been de- 
veloped for hdQCD which gives rise to the hope that hdQCD can be treated theoretically 
from the first principles. 

• We are very close to understanding of the parton density saturation 

In general, we think that the theoretical approach to hdQCD in a good shape now. 

3 HERA Puzzle: Where Are SC? 

The wide spread opinion is that HERA experimental data for > 1 GeV^ can be described 
quite well using only the DGLAP evolution equations, without any other ingredients such as 
shadowing corrections, higher twist contributions and so on ( see. for example, reviews [§] ). On 
the other hand, the most important HERA discovery is the fact that the density of gluons ( gluon 
structure function ) becomes large in HERA kinematic region |^ . The gluon densities extracted 
from HERA data is so large that parameter k ( see Eq.( |^) ) exceeds unity in substantial part of 
HERA kinematic region (see Fig. 2a ). Another way to see this is to plot the solution to Eq.( |^) ( 
see Fig.2b). It means that in large kinematic region k > 1 ( to the left from line k = 1 in Fig.2b), 
we expect that the SC should be large and important for a description of the experimental data. 
At first sight such expectations are in a clear contradiction with the experimental data. Certainly, 
this fact gave rise to the suspicions or even mistrusts that our theoretical approach to SC is not 
consistent. However, the revision and re-analysis of the SC , as has been discussed in the previous 
section, have been completed with the result, that k is responsible for the value of SC. 

Therefore, we face a puzzling question: where are SC?. Actually, this question includes, at least, 
two questions: (i) why SC are not needed to describe the HERA data on F2, and (ii) where are 
the experimental manifestation of strong SC. The answers for these two questions you will find in 
the next three sections but, in short, they are: SC give a weak change for F2 in HERA kinematic 
region, but they are strong for the gluon structure function . We hope to convince you that there 
are at least two indications in HERA data supporting a large value of SC to gluon density: 

1. xp - behaviour of the cross section of the diffractive dissociation (o"^^) in DIS; 

2. _ behaviour of F2 -slope ( ^g/^g'^ ). 
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Fig.2a 



Fig.2b 



Figure 2: Parameter k and new scale Q'q{x). 



4 SC for 



It is well known, that 7* - hadron interaction goes in two stages: (i) the transition from virtual 
photon to colour dipole and (ii) the interaction of the colour dipole with the target. To illustrate 
how SC work, we consider the Glauber - Mueller formula which describes the rescatterings of the 
colour dipole with the target pO[]: 



OTT-^ ^ ■' jQl J 

+ F^{xB,Ql). (8) 

where S{ht) = e "r^ is the target profile function in the impact parameter representation and 
I h{xb, Q'^) = crixB, r"^ = is the cross section of the dipole scattering in pQCD. 

One can see that Eq. ( ||) leads to 

F,ixB,Q') — ^,T.zjQ'R'. (9) 

However, we are sure that the kinematic region of HERA is far away from the asymptotic one. 
The practical calculations depend on three ingredients: the value of R^, the value of the initial 
virtuality Qq and the initial F2 at Qq. We fix them as follows: = 10 GeV~^ which corresponds 
to "soft" high energy phenomenology fill, Ql = lGeV^ and F2*"^"*(xb, Q^) = F^^^'^\xb,QI 



IGeV'^). Therefore, the result of calculation should be read as "SC for colour dipoles with the 
size smaller than r'j_ < A/lGeV"^ are equal to ..." 
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Figure 3: SC for xG{x,Q'^) and F2 in the HERA kinematic region 



From Fig. 3 one can see that the SC are rather small for F2 but they are strong and essential for 
the gluon structure function. It means that we have to look for the physics observables which 
will be more sensitive to the value of the gluon structure function than F2. 



5 Xp- dependence of a 



DD 



One of such observables is the cross section of the diffractive dissociation and, especially, the 
energy dependence of this cross section. 

Data: Both HI and ZEUS collaborations [§] found that 

« ^ ' (10) 
where Ap = ap(0) — 1 and the values of ap(0) are: 



HI g : ap(0) = 1.2003 ± 0.020(stat.) ± 0.013(sys) ; 



• ZEUS : «p(0) = 1.1270 ± 0.009(stat.) ± 0.012(sys) . 

It is clear that the Pomeron intercept (q;p(0)) for diffractive processes in DIS is higher than the 
intercept of "soft" Pomeron |]TT] . 

Why is it surprising and interesting? To answer this question we have to recall that the 
cross sections for diffractive production of quark-antiquark pair have the following form in pQCD 



dalnir - q + q) f"^ dkl Us xpG{xp,^^, 
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« L X ^ u ^; (11) 



dxpdt JQl kj^ kj_ 

fc2 



Wn-^ (^*^n^n) dk"^ [ as Xp G{xp, j^) 

. . oc / ^ X A ^ ^ / . (12) 

dxpdt JQl Q2 ^2 V ; 
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From Eqs.( [Tl|) and ( |T2D you can see that /cj_ integration looks quite differently for transverse 
and longitudinal polarized photon: the last one has a typical log integral over fc^ while the former 
has the integral which normally converges at small values of k^. We have the same property for 
production a more complex system than gg, for example qqG |jl5[. Therefore, we expect that the 



diffractive production should come from long distances where the "soft" Pomeron contributes. 
However, the experiment says a different thing, namely, that this production has a considerable 
contamination from short distances. How is it possible? As far as we know, there is the only one 
explanation: SC are so strong that xG{x, k\) oc k\R^ (see Eq.( ^ ) Substituting this asymptotic 
limit in Eq.( |lT]) one can see that the integral becomes convergent and it sits at the upper limit 
of integration which is equal to k\ = Ql{x). 
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Finally, we have 

--SS - (.PG(.P.QS(.P)))^ X ^ (13) 

The calculation for ^^^^I'jl^ is given in Fig. 4 for HERA kinematic region using Glauber-Mueller 
formula [|l^ for SC. Taking into account that Q'^{x) in Fig. 2b can be fitted as Qg(x) = I GeV^ [f^) 



with A = 0.54 and Xq = 10~^ we see from Eq.( |T^) and Fig.4 that we are able to reproduce the 
experimental value of ap(0) and conclude that the typical which are dominant in the integral 
is not small {k^ ~ 1 — 2GeV^ |]T3[ ). For Golec-Biernat and Wiisthoff approach, which we will 
discuss late, A = 0.288 and the value of typical k^ turns out to be higher. 





Figure 4: Effective Pomeron intercept for gluon structure function calculated using Glauber-Mueller 
formula for SC 
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6 The Q^- dependence of the F2 - slope. 



Data: The experimental data ||TB[ for F2 - slope dF2{x,Q'^) / dlnQ"^ are shown in Fig. 5a 
(Caldwell plot). These data give rise to a hope that the matching between "hard" ( short distance) 
and "soft" (long distance) processes occurs at sufficiently large since the F2-slope starts to 
deviate from the DGLAP predictions around ~ 5 — SGeV^^. 

F2-slope and SC: Our principle idea, as we have mentioned in the beginning of the talk, is 
that matching between "hard" and "soft" processes is provided by the hdQCD phase in the 
parton cascade or, in other words, due to strong SC. The asymptotic behaviour of F2 cx Q^F? 
for < Q2^a;) leads to dF2{x,Q'^) / d\nQ'^ oc Q'^R^ at < gs^^,) (gee Eq.( |) ) and this 
behaviour supports our point of view ||16| ITtI . 

However, we have two problems to solve before making any conclusion: (i) the experimental 
data are taken at different points (x, Q"^) and therefore could be interpreted as the change of 
x-behaviour rather than one; and (ii) the value of F2-slope is quite different from the value 
of F2 while for the asymptotic solution it should be the same. Therefore, we have to calculate 



F2 - slope to understand them. The result of calculation using the Glauber-Mueller formula |]16 
is presented in Fig. 5b. One can see that (i) the experimental data shows rather - behaviour 
then the x- dependence, which is not qualitatively influenced by SC; and (ii) SC are able to 
describe both the value and the Q^-behaviour of the experimental data. Fig. 5b shows also that 
the ALLM'97 parameterization [|18|, which can be viewed as the phenomenological description of 
the experimental data, has the same features as our calculation confirming the fact that the data 
show - dependence but not x-behaviour of the F2-slope. 



7 Golec-Biernat and Wiisthoff Approach 

Golec-Biernat and Wiisthoff suggested a phenomenological approach which takes into account 



the key idea of hdQCD, namely, the new scale of hardness in the parton cascade. They use for 
7*p cross section the following formula [|l^ 

(^tot{l*p) = I d^r^J^ dz\'^{Q^-r^,z)\^ (rtot{rl,x) ■ (14) 



a(x,r^) = ao { 1 - e ^ } ; (15) 



X 



-A 



R\x) = l/Qtix) with Qi{x) = — . (16) 

V Xo / 

Extracting parameters of their model from fitting of the experimental data, namely, (To=23.03 
mb, A = 0.288, Ql = IGeV"^ and Xq = 3.04 10""^ , they described quite well all data on total and 
diffractive cross sections in DIS (see Fig.6 ). 
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Figure 5: F2-slope: experimental data (Caldwell plot ) (Fig. 5a) and calculations using Glauber-Mueller 
formula (Fig. 5b) . 

8 Why Have We Only Indications? 



The answer is: because we have or can have an alternative explanation of each separate fact. 
For example, we can describe the F2 -slope behaviour changing the initial x-distribution for the 
DGLAP evolution equations pO|. Our difficulties in an interpretation of the experimental data is 



seen in Fig. 6a where the new scale Ql{x) is plotted. One can see that Ql{x) is almost constant in 
HERA kinematic region. It means that we can put the initial condition for the evolution equation 
at Qo =< Qoi^) > where < Ql{x) > is the average new scale in HERA region. Therefore, SC 
can be absorbed to large extent in the initial condition and the question, that can and should 
be asked, is how well motivated these conditions are. For example, I do not think that initial 
gluon distribution in MRST parameterization |20|, needed to describe the F2 - slope data, can be 
considered as a natural one. 



9 Summary 

We hope we convinced you that (i) hdQCD is in a good theoretical shape; (ii) the hdQCD region 
has been reached at HERA; (iii) HERA data do not contradict the strong SC effects; (iv) there are 
at least two indications on SC effects in HERA data: behaviour of F2 slope and xp behaviour 
of diffractive cross section in DIS; and (v) the HERA data and the hdQCD theory gave an impetus 
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ure 6: (ytot{l*p) amnd the ration a jotot for DIS in Golec-Biernat and Wiisthoff model. Vertical 
line in Fig.6a is the Q\{x) given by Eq.(16). 
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for a very successful phenomenology for matching "hard" and "soft" physics. 

We would like to finish this talk with rather long citation: "Small x Physics is still in its infancy. 
Its relations to heavy ion physics, mathematical physics and soft hadron physics along with a 
rich variety of possible signature makes it central for QCD studies over the next decade" ( A.H. 
Mueller, B. Miiller, G. Rebbi and W.H. Smith "Report of the DPF Long Range Planning WG 
on QCD"). Hopefully, we will learn more on low x physics by the next Ringberg Workshop. 
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